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Abstract 

A very simple model for cell swelling by osmosis is introduced, resulting in a parabolic free 
boundary problem. In case of radially symmetric initial conditions, it is shown that the model 
can be viewed as a gradient flow involving entropy, surface area and the Wasserstein metric. 
This observation is used to construct solutions and explain the presence and nature of osmosis. 
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1 Introduction 

In this article, a very simple model for osmotic cell swelling is introduced, and studied using a gra- 
dient flow approach. Observations of osmotic effects indicate that cell membranes are permeable to 
water, but impermeable to other molecules. Therefore, it is believed that water transport through cell 
membranes is facilitated by proteins acting as water channels, generally called 'aquaporins' ( HCMIO . 
Although there is growing evidence that aquaporins are proteins embedded in the cell membrane 
dMO. existence has been shown for only a few types of cells. The exact nature of aquaporins remains 
uncertain, although aquaporins are believed to be proteins (" HLVZIO . 

1.1 A model for osmotic cell swelling 

Quantitative modeling of water flux through membranes by aquaporins seems to have started with 
HLVZL In this paper, the water flux was assumed to be proportional to the difference in molar free 
energy inside and outside the cell. However, the model did not predict the eventual slowing of cell 
swelling. Several other models using this assumption have been considered since, but no model 
succeeded to provide evidence for either the permeability or aquaporin model when comparing 
predictions and observations. Furthermore, most models assumed spatially constant concentrations 
inside and outside the cell. 

Recently, Pickard (0), made a first attempt to develop a model for cell swelling, based on 
observations of Xenopus oocytes submitted to hypoosmotic shocks, that is, sudden exposure to low 
concentrations of solute. The proposed model uses proportionality of free energy differences and 
water flux as a starting point. In this model, the surface tension of the cell membrane is not taken 
into account. 
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Unfortunately, the justification of dropping surface tension does not seem to be valid if the cell is 
close to equilibrium. Although it should be noted that in most test tube situations the cell will burst 
long before it reaches equilibrium, the implications of assuming zero surface tension should be 
investigated. Furthermore, it might be insightful to compare a mathematical model with aquaporins 
with a similar model where the cell membrane has a constant permeability. Except for the last 
section, the constant permeability model will be considered. 



1.2 Mathematical formulation 

The model consists of two main ingredients: a time-varying domain with diffusion of mass inside. 

^ = kAu. (1.1) 
at 

The free boundary is non-permeable to the diffusing mass, which is basically a Neumann condition 
for a free boundary: 

- = uv. (1.2) 

on 

Here, u(x,t) is the concentration of solute, and v is the normal velocity of the cell membrane, where 
the normal n is pointing outward. It is comparable to the standard Neumann boundary condition 
since it requires that there is no flux of solute through the boundary of the free domain. This can be 
seen by introducing the additional variable v(x, t), the velocity of the diffusing mass. Since (11.11) 
implies that v = — k— , (11.2b is equivalent to 

uv-n = uv, (1-3) 

that is, normal velocity of the boundary and the normal component of the velocity of mass at the 
boundary should match. Note that this condition (11.2b does not depend on how the membrane reacts 
to u. 

In the present model, the normal velocity of the membrane is caused by the absorption of an 
incompressible solvent, usually water. Assuming constant permeability for the cell membrane, the 
amount of solvent absorbed through a portion of surface is proportional to the difference in the 
Helmholtz free energy A inside and outside the cell, which is given by 

AA = AP + Au, (1.4) 

where AP and Au denote the pressure difference and concentration difference accross the cell 
membrane, respectively. The pressure difference is due to the surface tension of the cell membrane, 
resulting in a term proportional to the mean curvature of the membrane. Since it is assumed that 
there is no solute outside the cell, the concentration difference is equal to just the value of u at the 
boundary. Since the solute is assumed to be incompressible, the normal velocity of the boundary 
should be proportional to the amount of solvent absorbed by the cell. Combined with the above 
model for the absorption of sovlent, this resulting equation for the normal velocity of the cell mem- 
brane is 

v = 1 H + fiu, (1.5) 

in which 7 and (3 are positive parameters. The sign convention for H is according to n, which means 
that H = — for a sphere of radius r in R™. Note that the first term will tend to make v smaller, 
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whereas the second term tends to make v bigger. This is in accordance with the intuition of the 
model: the surface tension of the cell membrane will tend to make the cell smaller, whereas the 
osmotic force tends to make the cell bigger. 

A number of parameters from this problem can be eliminated by changing variables. Looking at 
the units of u, x and t, it follows that ac and (5 have the same dimension. Therefore, at most one of 
the two can be set equal to 1 in general by changing units. 

Since the total mass is a conserved quantity, it can also be regarded as a parameter, say i?. 
Changing units for x first, the total mass can can be set equal to 1. Next, changing units for u 
and x simultaniously, (3 can be made equal to 7 without changing the total mass. Finally, scaling t, 
7 = j3 = 1 can be realized. This results in new, scaled variables 



R\ 5TT „_j / R \ H+T 



(1.6) 



It is easy to check that w and p satisfy (11.11) . (|1.2I) . (11.51) . with coefficients 7 = 1 and (3 = 1, and 
f w d£ = 1 . Note that the change of coordinates also implies a change in the coefficient ac. 
It follows that one can restrict attention to the problem 

du 

— = kAu, for x e Q(t), t > 0, 

—k-^— = uv, for x G dQ(t), t > 0, (L7 ' ) 

an 

v = H + u, for x E dCl(t), t > 0, 

with the integral of u equal to 1 . The remaining parameter n, which cannot be scaled out, can now 
be regarded as a measure of how fast the dynamics is inside the domain, compared to the movement 
of the boundary. 

Looking at (11.71) . a few observations can be made. First of all, the diffusion equation can be 
rewritten as a system of two equations with two unknowns: 




-ac , 



dt 



Vu 

u 

— div(iiv) 



Since v is regarded as the velocity of mass particles at a certain position and time, the first equation 
is, up to multiplication by u, Fick's law of diffusion, which can be regarded as a modelling law 
describing how mass moves in reaction to the spatial variations of the concentration u. In constrast, 
the second equation is conservation of mass. The boundary condition, which was given in terms of v 
in (11.31) . reduces to the requirement that the normal part of the velocity of the mass and the velocity 
of the boundary should math whenever u is positive. From a modeling point of view this formulation 
as a system is more natural: each equation and each unknown has a physical interpretation. In the 
classic approach, one of course eliminates the velocity and writes the diffusion equation with only 
u in it. This representation as a system representation can also be found, for instance, in jQl §2]. 
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Here, the Porous Media Equation is reformulated in the same way, by replacing only Fick's law by 
another law. 

A similar observation can be made for the time-varying domain. The splitting between a mod- 
eling equation and a continuity equation has already been done. The modeling equation is, in this 
case, the expression for the normal velocity v. The continuity part is not written explicitly, but in 
words, it is the remark that v is the normal velocity of dCl(t). This can be put into a formula in a 
number of ways, depending on the smoothness of dQ(i). 

Summing up, the problem (11.71) can be viewed as a system of four relations, each with a physical 
interpretation: 

• Fick's law of diffusion, giving an expression for the velocity of particles, 

• Conservation of mass, relating the change in u to the velocity, 

• The modelling law v = H + u, 

• v being the normal velocity of <9f2. 

As explained above, these four relations can be divided in two groups in two different ways: the first 
two describe the evolution of u, whereas the latter describe Vt. On the other hand, one can argue that 
the first and third contain the actual model that is being studied, whereas the second and fourth are 
merely equations relating the evolution of u and Vt to the functions v and v, respectively. These four 
equations will play a key role in what follows, as will the two subdivisions made here, which are 
natural from a modeling point of view. 

Note that the problem is almost a superposition of two well-known problems: the diffusion 
equation and the mean curvature flow. However, there is a nonlinear coupling, which consists of 
two parts. First of all, u is required to be supported inside Vt, and secondly an osmotic boundary 
term appears in the evolution of Vt. It will turn out that there is a close relation between these two 
couplings; in some sense, the latter follows from the former, as will be shown later. 

For both the diffusion equation and the mean curvature flow, variational approaches have been 
developed. A well-known paper about a variational approach to the Fokker-Planck equation, and 
the diffusion equation in particular, is HJKOII . It has turned out that this is a fundamental example 
of a flow that can be seen as a gradient flow in a metric space, for which a fairly general theory 
is developed in HAGSL with gradient flows in the space of probability measures, and diffusion 
in particular as an example. Beside the diffusion equation, many parabolic equations have been 
modeled as gradient flows in the space of probability measures with the Wasserstein metric, see for 
instance (D]|, which also treats the Porous Medium Equation, or [AJ. 

A variational approach for the mean curvature flow can be found, for example, in [ATWl or ULSL 
Although there is, a far as the author is aware, no theorem stating that the mean curvature flow is in 
fact a gradient flow in a metric space, the similarities between the discretization used in [ATWl and 
ULSI on the one hand, and the discretization used for the diffusion equation in [JKOJ on the other 
hand, strongly suggest that the mean curvature flow is a gradient flow. The functional playing the 
role of the metric in HATWI is, however, not a metric. Recalling the first variation of area formula, 
an appropriate metric would be related to the L 2 norm of the normal velocity, as is noted in HGKI . 
Unfortunately, this yields a trivial metric in general. 
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One way to avoid this problem and obtain a gradient flow with respect to an honest metric, is 
to restrict attention to the case of radial symmetry. Although this essentially reduces the boundary 
movement to a one-dimensional problem, the resulting problem is still rich enough to see how 
exactly the osmotic term in the normal velocity v arises from the restriction on the support of u. 

Another possible approach to the problem is to combine the time discretization techniques in 
ULSI and [JKOJ without making an explicit connection to gradient flows. Obviously, this is techni- 
cally much more involved. Although some of the results are easily generalized, taking the limit for 
the discretization parameter to zero is much harder than in either of the two individual cases. 

In the following sections, it will be shown that the radially symmetric version of (11.71) is a gradi- 
ent flow in a metric space, and solutions will be constructed using this observation. More precisely, 
in Section [2l a metric space and functional will be constructed for radially symmetric initial con- 
ditions, and fitted into the framework of [AGS], enabling the use of a theorem that guarantees the 
convergence of an abstract version of the approach used in [JKO], ULSH and HATWH . In Section [3l 
some finer properties of the space and functional from Section [2] will be studied, resulting in a theo- 
rem concerning, among other things, regularity of the curve constructed in Section [21 Convexity of 
certain functionals will be a special point of interest. Section |4]is devoted to the study of the differ- 
ential properties of the functional constructed in Section[2l making the conclusions of Theorem |3.9l 
more concrete. In Section[5]it is shown, using theory from Sections [3] and |4] that the construction in 
Section [2] indeed gives a weak solution of (11.71) . Finally, in Section|6l a model for aquaporins will be 
presented. It will be shown that a very small modification of the metric space is sufficient to obtain 
a gradient flow formulation for the resulting problem. 



2 Construction of a generalized minimizing movement 

In order to obtain a gradient flow solution for (|1.7I) under the assumption of radial symmetry, a 
metric space will be constructed to represent the domain and mass profile. It will be shown that this 
space, equipped with a weak topology, and a suitable functional fit the theory in HAGSl §2]. More 
precisely, a so-called Generalized Minimizing Movement or GMM will be constructed. The idea of 
this concept is a (formal) generalization of the Euler backward approximation of a gradient flow in 
R n . 

The gradient flow of a function <p on W 1 can be constructed by an Euler backward scheme. Time 
is discretized by setting tk = kh for some small h > 0, and the equation 

V0(x fc+1 ) + Xk+1 ~ Xk = 0, (2.1) 

is solved for every timestep. Setting x(kh) = x k , this is an approximation of the gradient flow of (p. 
Note that (12.11) is the optimality condition for the minization problem 

This last problem can easily be generalized to the setting of a metric space by replacing \y — x\ 2 
with d 2 (x, y). It is expected that this will, in some sense, give an approximation of a gradient flow. 
In general, (12.21) might not have a unique solution, or no solution of all. Moreover, one must ask 
whether the limit for h 4- exists, in what sense, and whether this limit is unique. 
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The concept of a (Generalized) Minimizing Movement is defined in HAGSl Definition 2.0.6]. 
The starting point is a partition of the time interval (0, +00) 

V T := {0 = t° < 4 < ■ ■ • < t k T < ■ ■ ■ } (2.3) 

associated to a sequence r = {r fc } feeN of positive time steps with 

00 

lim t k = J2 T i = +°°- ( 2 - 4 ) 
i 1 

and |r I := sup fc {rfc} < +00. Additionally, the time intervals 

I k T :={t k T -\t k T ] (2.5) 

are introduced. 

Given a complete metric space (S, d) and a functional : S — > (—00, +00], a discrete solution 
associated to the partition V T is a map {7 T : (0, +00) — > S such that U T is constant on I k for every 
G N, and U T (t k ) is a minimizer of the functional 

V H- ^(CTrCt*- 1 ), V) := 0(V) + J-d 2 ^^- 1 ),^). (2.6) 

2r fc 

A curve [0, +00) — > S is called a minimizing movement for starting at « G 5 if for every 
partition r with |r | sufficiently small there exists a discrete solution [7 T such that 

lim (j>(y T (0)) = (j)(u ), limsupd(!7 T (0), w ) < +00, (2.7) 

l T U-° |r|4.0 

and U T (t) — >■ w(t) for alH > as |r| J, 0. In this definition, convergence of the discrete solutions 
with respect to the metric may be replaced by convergence with respect to a weaker topology. 
This allows for the use of a topology that may enjoy better compactness properties than the metric 
topology. 

A generalized minimizing movement is defined analogous to a minimizing movement, but with 
(|2.7I) required only along a sequence Tj of parititions. That is, a curve [0, +00) — > S is a generalized 
minimizing movement for starting at uq E S if there exists a sequence of partitions Tj with 
lirn^oo I I = and discrete solutions U Tj such that 

lim 0(I7 T .(O)) = (f>(uo), limsupd(I7 T .(0),Mo) < +00, (2.8) 

and U Tj (t) — > u(t) for all t > as j — >■ 00. Again, metric convergence may be replaced with 
convergence with respect to a weaker topology. 

Examples of the construction of a GMM can be found in HJKOH . HATWI and HLS1 . where this 
approach is applied. In this section, a metric space (X n , g) and a functional $ on X n will be intro- 
duced. Moreover, it will be shown that GMM's exist for initial conditions that have finite $-value 
by applying HAGS I Proposition 2.2.3], which is proven in §3 of the same book. This results in the 
following theorem. 
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2.1 The space of balls 



By radial symmetry, the free domain will be a ball centered at the origin, which means it can be 
represented by a non-negative real number: the space of balls centered at zero is C n := [0, +00). 
The abbreviation B r for B r (0) will be used throughout the paper. Following HATWL ULSL [|Bl and 
many others, in order to obtain the mean curvature term the perimeter P n (r) = nu n r n ~ l of dB r 
will be used. 

Although metrics on C n are trivial to construct and study, it is interesting to view them as the 
restriction of metrics on a suitable space of domains to the space of balls. 

A natural way to construct such a metric is described in HGKl §2] . One could consider all families 
{E t : t 6 [0, 1]} of domains with smooth boundary, where E and E\ are prescribed, such that 

U dEt x {t} (2.9) 

te[o,i] 

is a smooth, n-dimensional manifold. If v t is the normal velocity of dE t , the integral 



\JdEt 



2 



v 2 t m n - 1 ) dt, (2.10) 



is positive whenever E ^ E\. One would like to define the distance between E and E\ to be the 
infimum of (12.101) over all such families {E t : t 6 [0,1]}. Unfortunately, this infimum is in general 
not positive, even if E and E\ are different. However, if the E t are restricted to be balls centered at 
zero, the infimum can be computed explicitly: 



d(r ,ri) := / V ' P n (p) dp 

J (Vn,n ) 



'(ro.ri) n + 1 



r i 2 ~ r o 2 



n + l 

(2.11) 



where E = B ro and E(l) = B ri . The fact that the above infimum is positive when restricted to C n 
is the main reason to consider only the radially symmetric case. 

In the general situation, the integral of a distance over the symmetric difference of two domains 
in llATWl §2.6] and PKOl (1)] instead of d: 

d 2 (B ro , B ri ) = 2 / d(x, dB ro ) dx . 

J B ro AB ri 

However, this is not a metric: it is not symmetric and does not satisfy the triangle inequality. 

If it does not matter which one is used when constructing solutions, one would expect that, at an 
infinitesimal level, the two are the same. This is indeed the case: by evaluating d for balls, 

r / r n+l r n+l r ry \ 

d 2 (r , n) = 2 / d(x, dB ro ) dx = 2nu n + - — ) , (2.12) 



B ro AB ri 



n+l n(n + 1) n 



so that, 

-77<l{ro,rt 



dt 



t=o 




r'^/nwnrr 1 ^ I / \v\ 2 dU^ 1 ) (2.13) 
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since v = r' for balls. Intuitively, this means that the two metrics have 'the same infinitesimal struc- 
ture' . What this means exactly, and what the implications for gradient flows are is being investigated 
by the author. Here, (12.1 II) will be used as the metric on C n . 
Clearly, the map 

1\/nUJ n n+l 

i n : C n -> 0, oo : r h- -^-^ r — (2.14) 

n+l 

is an isometry, which will be used in computations. One can see immediately that (C n , d) is a 
complete metric space. Moreover, (|2.1 II) generates the standard topology on C n , which means in 
particular that P n is continuous with respect to d. Finally, — > r is equivalent to £ n (B rk AB r ) — > 
as k — > oo. 



2.2 Mass profiles 

The mass profile will be a radially symmetric nonnegative function with a fixed integral. Identifying 
a profile with a function u on (0, oo), fixing the total mass is equivalent to requiring 

u(p)P n (p)dp= 1 (2.15) 

Hence, an appropriate space for the solute profile is the space L n of radially symmetric proba- 
bility density functions on R™, which is a metric subspace of L 1 (M ri ). Alternatively, it can be seen 
as a subspace of ^((0, oo), P n ), the space of integrable functions on (0, oo) with weight P n . 

The topology on this space will be the weak topology inherited from L 1 (P n £ 1 ), by definition, 
this means that u k — u if and only if 

/*oo roo 

/ u k (p)g(p)P n (p)dp^ u(p)g(p)P n (p)dp Vf/GL oo ((0,oo)). (2.16) 
Jo Jo 

The following compactness result is an easy consequence of the Dunford-Pettis theorem jAFPl 
Theorem 1.38]. 

Theorem 2.1 (Relative sequential compactness). Let {w/c}fceN be a sequence in L\ such that 

• for any e > 0, there exists an R > such that u k (p)P n (p) dp < e for all k E N, 

• there exists a nondecreasing function g : [0, oo) — > [0, oo) with superlinear growth and a 
constant C E R such that 



/■oo 

/ g(u k )P n ML x <C. (2.17) 
Jo 



for all jfe e N. 

Then {M fc } fceN has a subsequence {uk(j)}jen converging weakly to u E L n . 

Proof. By ||AFPI Proposition 1.27], {wfejfcgN is equi-integrable. Therefore, using the Dunford-Pettis 
theorem jAFPl Theorem 1.38], it has a subsequence {uk^jem that converges weakly to some u E 
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L 1 ((0, oo), P n ). Since u k (p) > 0, it follows that u(p) > for /^-almost all p > 0. Moreover, since 
the constant function is certainly in L°°((0, oo)), 



POO POO 

/ u(p)P n (p) dp = lim / u k (p)P n (p)dp 
Jo k ^°°Jo 



1. 



which means that u 6 L\ . □ 

Note that all that needed to be shown is that L\ is closed as a subset of -^ 1 ((0, oo), P n ). 

The metric on L\ set will be the Wasserstein distance, which is based on the optimal trans- 
portation problem. Given two mass profiles u and w, one can consider all measurable maps T : 
[0, oo) — > [0, oo) such that the profile u turns into w if all mass is transported according to the map 
T. Mathematically, this means that the push-forward T#u of u under T, defined by 

/ u(p)P n (p)dp= ! T # u(p)P n (p)dp (2.18) 

J A JT{A) 

for every Borel set A C [0, oo), must be equal to v. The Wasserstein metric is defined as 



W 2 (u,v):=( inf / \x-T(x)\ 2 dx) , (2.19) 

\T:T # u=v J Q J 

which can be interpreted to be the minimum amount of work needed to transport mass from a 
profile u to another profile w. For an extensive introduction of the optimal transportation problem 
and the Wasserstein metric, see HAGSl §6-7]. The heat equation as a gradient flow with respect to 
the Wasserstein metric is treated as an example in this book. 

The optimal transport problem in this setting is essentially the same as the one-dimensional 
optimal transportation problem. Therefore, the result of [AGS, Theorem 6.0.2] can be carried over 
to this setting. Writing 

F u {p) := [ uPndC 1 (2.20) 
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the optimal transport map t™ : [0, oo) — > [0, oo) and Wasserstein metric are given by 

t: = F~ l oF u (2.21) 

poo i* 

W*(u,w)= \p-F w 1 oF u (p)\ 2 u(p)Pn( P )dp= iF-^-F-^fda (2.22) 
</o Ao,i) 

whenever F u and F w are invertible. If F u or F w is not invertible, the same identities hold if the 
inverse is replaced by the pseudo-inverse 

F-\a) := sup {p > : F u (p) < a} , a G [0, 1]. (2.23) 

This observation will be useful later for studying convexity properties of the Wasserstein metric. 

In HAGSl Proposition 7.1.3], it was shown that the Wasserstein metric is lower semicontinuous 
with respect to the topology of narrow convergence. Since the weak L 1 topology is stronger, it 
follows immediately that W2 is also lower semicontinuous with respect to this topology. 

Note that L\, equipped with Wasserstein metric is not a complete metric space, a sequence 
of profiles concentrating at the origin has no limit in L\. Moreover, the weak L 1 topology is not 
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necessarily weaker than the topology generated by the Wasserstein distance. It will turn out that this 
will not be a problem due to the choice of the functional, which will be defined below. 

As in HJKOH . the entropy functional will be used. Given a mass profile, its entropy is defined as 

poo 

/ u(p) log u(p)P n (p) dp, 
Jo 

which is, up to a constant, Boltzmann's Entropy. Since the exact formula of the integrand is of little 
importance to the analysis, z log z will be replaced by a more general function. 

Definition 2.2 (Internal energy). Let / : [0, oo) — > R be a function such that 

• /(0) = 0, 

• lim^o f(z) = 

• / is strictly convex, 

• lim^ +00 ttf- = +oo, 

• / is continuously differentiable on (0, oo). 
Then 

POO 

F{u):= / f(u(p))P n (p)dp (2.24) 
Jo 

is called the internal energy of u. 

Remark 2.3. By continuity of / and superlinear growth, it follows that / has a global minimum, 

say 

/(*) > -/o. (2.25) 
Note that the minimum of / does not need to be negative, / may be increasing. 



Since the integrand is assumed to be convex, the internal energy defines a weakly lower semicon- 
tinuous functional on L\. This has been proven in a much more general setting in ||AFP1 Theorem 
2.34] by writing / as the supremum of countably many linear functions. 

In what follows, the following auxiliary function will be used: 

/: [0,oo)^ [0,oo);^ <f W - /W - " 2> °' (2.26) 

0, if z = 0. 

Note that, since / is convex, / is increasing, nonnegative and continuous. If f(z) = z log z, f(z) = 
z. 
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2.3 The variational formulation 



The two spaces from the previous subsections can be combined into one space. This requires some 
caution, because the mass is supposed to stay inside the free domain. Therefore, the following space 
will be used. 

X n :=^(r,u)eC n xL 1 n : uP n dC x = ^ (2.27) 

Note that the support of u need not be the whole ball B r . 
The metric on this space will be 

g((r, u), (s, w)) := (d 2 (r, a) + * . (2 .28) 

Moreover, since X n is a subspace of a product, it inherits a topology from C n and L\: weak con- 
vergence in X n , denoted by (r fc , u k ) — 1 (r, w) is equivalent to r k — > r and u k — 51 it, where u k ^ u 
denotes weak convergence, characterized by (12.161) . 

An obvious question is if the compactness properties of the spaces C n and L\ carry over to X n . 
Clearly, using Theorem 12. II and the Heine-Borel theorem for M 1 , suitable sequences C n x L\ have 
convergent subsequences. The question is, however, if the condition that all mass should stay inside 
the varying domain can be carried over from a sequence in X n to its limit, if it exists. The following 
theorem states that this can be done. 

Theorem 2.4. X n is closed as a subspace of C n x L\ in the weak topology. 

Proof. Suppose that (r k , u k ) — ^ (r, u) in C n x L\ with (r fc , G for all k E N. By the Dunford- 
Pettis theorem HAFPl Theorem 1.38], {nfc}fc e N is equi-integrable. Let then e > 0, and pick 6 > be 
such that 

£«(A)<5^ [u k (x)dx<e (2.29) 



for all k E N and all Borel sets A C R n . Choose also N E N such that, w„(r£ - r n ) < <5 for all 
k > N. Then, 

u k (p)P n (p) dp< u k (p)P n (p) dp + / u k (p)P n (p) dp < e, (2.30) 

Jmax{rfc ,r} •/ r 

because the first integral is zero by definition of X n and the second integral is smaller that e by 
(12.291) for all k> N. Taking limits for k -> oo, it follows that 



w(p)P n (p) dp < e, (2.31) 
for any e > 0. Then the integral must be zero, which means that (r, u) E X n . □ 

The functional that will be used is a combination of the perimeter and the internal energy: 

$(r,«) :=P n (r)+F(u). (2.32) 
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In order to obtain coercivity, an assumption on the integrand / in relation to the dimension is made: 

lim z~»/(z) = 0. (2.33) 

This condition can be interpreted as a coercivity condition. Without it, very large balls with a con- 
stant mass profile have very low $-values. This would, in the end, result in the domain growing 
indefinitely large. A nice example is the case where the dimension n is 1 and f(z) = z log z. Note 
that, in particular, the choice f(z) = zlogz does satisfy the assumption, provided that n > 1. In 
order to ensure coercivity, it will be assumed that n > 1 and (12.331 ) holds in the remainder. 

It turns out that, once the assumption is made, the sublevels of $ are complete with respect to q 
and the weak topology is the same as the metric topology on sublevels. Before proving this, some 
basic properties of $ will be shown. 

Lemma 2.5. Given r G C n , h 4 $(r, u) takes its unique minimal value 

$ Q (r) := nu n r n - 1 + u n r n f ( (2.34) 



at u = — ^ Ym rl- The function $ G is coercive, in the sense that 

• lim r ^ +00 $ (r) = +oo, 

• lim r;o $ G (r) = +oo, 

• t i — y $ (r) has a unique global minimum. 

Proof. Since / is convex, Jensen's inequality flFj Theorem 2.4.19] implies that u = — ^X(o,r) 
minimizes u) for given r > 0. Superlinear growth of / implies that lim r | ^o(^) = +oo. 
Moreover, lim,..^ $ D (r) = +oo using (12.331) . Since 

is strictly increasing, the second term in (12.341) is convex. Since the first term is obviously also 
convex, it follows that r i-)- $ (r) has a unique global minimum. □ 

Using this coercivity lemma, the following two properties of sublevels of $ can be shown. 

Theorem 2.6. For any M E R, $ _1 (-oo,M] := {{r,u) e X n : $(r,w) < M} is weakly 
sequentially compact. 

Proof. Let {(r^, Uk)}keN be a sequence in oo, M). By Lemma [231 {r^} is uniformly bounded 

away from and +oo, say c < < C, which means it must have a convergent subsequence, say 
rk — > r > 0. Taking the corresponding sequence in Uk, it follows that 

poo 

/ max{/(u(p)), 0}P„(p) dp < M - ^C™- 1 + o; n C n /o (2.36) 

along this subsequence. Since also ut(p) = for any p > C, Theorem 12. II with g = max{/, 0} 
implies that there must be a further subsequence converging weakly to some u 6 L\. Taking 
the corresponding subsequence for r k , this results in a subsequence {(rw), Wfcfj))}jeN such that 
(r fc (j), Mfc(j)) — > (r, u) for some (r, u) 6 C„ x By 12.41 (r, w) G <Y n . Moreover, by continuity of 
P n and weak lower semicontinuity of J 7 , $(r, u) < lim inf $(r fc , < M. □ 
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Lemma 2.7. Let {(r k , u k )} k& ^ be a sequence in $ 1 (— oo, M]. Then (r fe , u k ) — ^ (r, u) if and only 
if (rjt, Wfc) — >■ (r, m) with respect to g. 

Proof. Suppose first that (r k , u k ) — 51 (r, «). Then, as above, «fc(p) = for p > C. Therefore, the 
corresponding measures on R n have uniformly integrable second moments. By definition, weak 
convergence in L 1 of u k implies narrow convergence of the measures p k to p. Applying HAGS I 
Proposition 7.1.5], it follows that u k converges to u with respect to the Wasserstein metric. Then 
(Vfc, Uk) — > (r, u) with respect to g. 

Conversely, assume (r k , — > (r, m). Then boundedness of r k , together with Theorem [2J] im- 
plies that (r k (j), u k (j)) must have weak limit point (r', u'). Clearly, this implies r' = r. By definition 
of weak* convergence of measures and weak convergence in L 1 , the measures p k] corresponding to 
u kj must weak* -converge to jj,', the measure corresponding to u'. On the other hand, since u k — > u, 
jj kj — ^ fi, the measure corresponding to u as well. This implies y! = //, which means that v! = u 
almost everywhere. □ 

The above analysis shows that the metric space 

X n = | (r, u) £ C n x L\ : «P n d/: 1 = 1 , (2.37) 

, 2 /„ „n , wi(u,w)y 



g((r, u), (s, w))= [d\r, s) + ^ y j (2.38) 
and the functional 

/>oo 

$(r, «) = P n (r) + = w"- 1 + / f(u(p))P n (p) dp (2.39) 

satisfy the requirements of [AGS* §2.1] if the problem of finding a GMM can be restricted to a 
sublevel of $. Lower semicontinuity of g and $, coercivity of $ hold on the whole of X n . Since 
on sublevels of $, the weak and metric topology coincide, strong compactness follows from weak 
sequential compactness. In particular, sublevels of $ are complete. 

By requiring that $ is finite for the initial condition, it is no problem to restrict the construction 
to a sublevel of $. Therefore, HAGS I Proposition 2.2.3] applies, which finishes the proof of 

Theorem 2.8. Let (r ,u ) E X n , such that $(r ,M ) is finite. Then a GMM for $ starting from 
(r , Mq) exists. 



3 Interpolation, Convexity and Uniqueness 

In this section, the convexity properties of the problem are studied. The goal is to prove that HAGS I 
Theorems 2.4.15 and 4.0.4] can be applied. This implies a number of properties of the GMM con- 
structed in the previous section, and eliminates the finiteness condition from Theorem 12.81 

The main assumption is convexity of the minimization problem solved in every time step of the 
discretized problem. 
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Assumption 3.1. For every w, v , V\ in (X, d), there exists an interpolating curve 7 : [0, 1] — > X 
such that 7(0) = v , 7(1) = v±, and the map 

v m- <f>(v) + ^<P(w,v) 
is (x + A) -convex along 7 for every < h < 3—. 

Remark 3.2. This assumption implies that for /i < , any two points are connected by a curve 
along which (13.11) is strictly convex. In particular, this means that minimizers are unique. In turn, 
this implies that for \r\ < 3— , discrete solutions associated to r are unique. 

The definition of A-convexity along a curve will be recalled below. As is the case with ordinary 
convexity, this assumption is only useful if v and V\ have finite 0-value, but w, vq and v\ can be 
restricted to an even smaller set, see HAGSl Assumption 4.0.1] for the details. 

Before showing that (X n , g) and $ satisfy this assumption, some finer properties of (X n , g) will 
be analyzed further. 



3.1 Absolute continuity and the metric derivative 

In order to study interpolating curves, and geodesies of g in particular, it is convenient to first study 
the notion of absolute continuity and metric derivative, as introduced in HAGSl §1.1]. Summarizing, 
a curve 7 parametrized on [0, 1] is absolutely continuous with respect to a metric d if there exists a 
function g G L 1 ([0, 1]) such that 

d(7(*),7(*))< / 9(r)dr (3.2) 

J s 

for all < s < t < 1. If 7 is absolutely continuous, the limit 

called the metric derivative of 7, exists for almost every r. Moreover, it is the smallest g that satisfies 
(13.21) . Note that existence of the limit does not guarantee absolute continuity of the curve. In this 
section, the notion of absolute continuity will first be analyzed for d and W 2 first. 

First, absolute continuity in C n is studied. Using the isometry (12.141) , it is easy to connect absolute 
continuity with respect to d to absolute continuity with respect to the Euclidean metric. If r h-> r(r) 
is absolutely continuous with respect to d, 

d(r(s),r(t)) < J \r'\ d (r)dr, (3.4) 

where |r'|d is the metric derivative of the curve with respect to d. By applying the isometry i n , the 
curve r 1 — v p(r) := L n {r{r)) is absolutely continuous in [0, 00). It follows that p(r) is differentiable 
almost everywhere, and that the metric derivative with respect to the Euclidean metric is just the 
absolute value | | of the ordinary derivative for almost every r. Since t n is an isometry, this 
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absolute value must also be equal to |r'|(r) for almost every r. An easy computation now shows 
that r(r) must be differentiable, and 



\r'Ur) = VPM^T) 



dr(r) 



(3.5) 



dr 

It follows that absolute continuity of r is absolute continuity of r as an IR-valued function, together 
with integrability of ^ P n {r{j))r' (r). 

There is also a useful characterization of absolute continuity involving a weak formulation. 
Clearly, if r t is smooth, 

£ 1>(j>,T)P n (p) dp = ^Ap n (p) d p+^)^(r(r),r))P„(r(r)). (3.6) 

for any smooth, radially symmetric test function ip. As it turns out, this equation is sufficient to 
characterize absolute continuity of r(r), and provides an expression for the metric derivative. 

Lemma 3.3. Let r ^ r(r) be a curve in C n . Then r(r) is absolutely continuous if and only if it is 
continuous and there exists a function g such that g{r)yj P n {r{r)) G L[ oc ((0, oo)), and 

00 r (r) dil)(o r) r 00 

/ ^^P n (p)dpdr = - / ff(T)^(r(r),r)P B (r(r))dr (3.7) 
o Jo dT Jo 



for all -0 G C£°([0,oo) x (0, oo)). In this case, |r'| d (r) = \g(r)\y/P n (r{tr)) for almost every t > 0. 

Proof. Suppose first that r(r) is absolutely continuous. Then |r'| d (r) = \^p-\\/P n ( r ( r )) G 
^ oc ((0, oo)) and, for any G C c °°([0, oo) x (0, oo)), 



Jo Jo dr ho J J h 

lim / — 

hio Jo h 



r°° 1 f r r ( T ) 
lim / - / ^(p,r)P n (p)dp 



(3.8) 

r(r+h) \ 

ifj{p,r)P n {p) dp dr 
Jo ) 

TT ^P-i>(r(r),r)P n (r(r))dr, 
dr 

if h > is small enough. 

For the opposite implication, assume that r(r) is continuous and g satisfies the equation. One 
would like to substitute ip(p, r) = X[t ,t 1 ]'P(p) w ^ tn 4> e C£°(IP> °°)) mto <I3.7I) , but this function 
is discontinuous. Let then C fc : (0, oo) — > [0, 1] be a sequence of smooth functions such that ( k — > 
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X[t ,ti] pointwise. Substituting 0(p)Cfc("7") for ip in (13.71) . one obtains 



to 



k— >oo 



r(r))^(r)P„(r(r)) dr = lim / g{r)<P{r{r))Ur)Pn{r{r)) dr 



roc r r \ T ) 

lim / Cfc(r) / 0(p)P„(p)dpdr 
k ^°°Jo Jo 

f oo , f f r(r) \ ( 3 - 9 ) 

hm^ C fe (r)-My o 0(p)P„(p)dpjdr 



r(ti) 



0(p)P n (p)dp 



K*o) 



if G C£°([0, oo)) for all to < £i- Using another approximation argument, this also holds for 
0(p) = P n (p)~2, which implies that 



d(r(* ),r(*i)) 



r(i x ) 



V^pJdp 



< 



/(r)v/P n (r(r))dr 



to 



(3.10) 



Since t and ^ are arbitrary, this implies that r(r) is absolutely continuous with metric derivative 

\r'\ d (r) = \9(r)WPn(r(r)). ' □ 

Next, absolute continuity in L\ is considered. For this, it is helpful to remember that L\ can be 
considered as a subspace of the space of probability measures on W" 1 . A study of absolute continuity 
of curves in the space of measures can be found in HAGS I §8.3]. The main theorem HAGS I Theorem 
8.3. 1] from this section states that, essentially, the metric derivative of a curve of measures p(r) 
can be found by solving the continuity equation 

dp{r) 



dr 



+ div(v(r)p(r)) = 



(3.11) 



in distributional sense. It is shown that, if r i-> p(r) is absolutely continuous, the metric derivative 
at any time to is the minimal L 2 (p(r))-norm of solutions v(r) of (13.1 II) . Note that v(r) can be 
thought of as the velocity of the mass at a certain position and time. In particular, \\v(r) H^^-j) can 
be thought of as the kinetic energy. 

In order to specialize [AGS, Theorem 8.3.1] to L\, a particular space of test functions is needed: 
the space of functions ijj : [0, oo) x (0, oo) such that <p : R n x (0, oo) defined by </>(x, r) = ^(M? r ) 
is C°° and compactly supported will be denoted by C™([0, oo), (0, oo)). 

Theorem 3.4. Let r h-> u(r) be a weakly continuous curve in L\. If r ^ u{r) is absolutely 
continuous, there exists aBorel function v : (p, r) !->■ v T (p) such that v T G L 2 ((0, oo), u(r)P n ) with 
\ML^ao,oo),u(r)Pn) < KI( T ) for almost every r and 



oo poo 



V dr 



djKf), r) 
dp 



v T (p) u(p, r)P n {p) dpdr = 



(3.12) 



for all ijj G C%° r ([0, oo), (0, oo)). Conversely, if there exists v with v T G L 2 ((0, oo), u(r)P n ) such 
that (13.121 ) holds for every ip G C^.([0, oo), (0, oo)), then r i— > u(r) is absolutely continuous and 

W\(t) < \\v T \\ L 2 moo)Mr) p n) for almost every r > 0. 
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Proof. Suppose first that r i-)- u(r) is absolutely continuous. By [AGS, Theorem 8.3.1], there exists 
a Borel vector field v : (x,r) >->■ v T (x) such that v r G L 2 (]R n ,ii) with ||v T ||x,2QRn )U ) < |w'|(r) 
for almost every r G (0, oo). It follows that for such a vector field, HAGSl (8.3.8)] holds with 
restricted to radially symmetric, smooth, compactly supported functions on W 1 x (0, oo). Under this 
restriction, HAGS I (8.3.8)] becomes invariant under rotations, that is, if v T (x) = v(Rx) with R a 
rotation in W 1 , v also satisfies HAGSl (8.3.8)] for radially symmetric 4>. Hence, it may be assumed 
that v r is radially symmetric for almost every r > 0. Let ip G C~ ([0, oo), (0, oo)) be given, and set 
<f>(x, r) := ip(\x\,T). Then, using llAGSl (8.3.8)], 

° = r L (^T^ + ^tW^^x, t))) u(x) dx dr 

f dtp(x ,r) ^ / v f x \ JL\ ^'(IxIt) ) u(x) dxdr 

in V OT \ \X\ I ) 







(3.13) 



oo /-oo 



JO 



+ (v T (pei),ei)^'(p,r) u{p)P n {p) dpdr, 



V Or 

which means that v T (p) := (v T (pei), ei) solves (13.121) . Moreover, by construction, 

lkr|U 2 ((0,oo),«(r)P„) < ||v t ||l2(IR«, u ) < \u'(t)\ (3.14) 

for almost every r > 0. 

For the converse implication, assume that v T satisfies (13.121) for every ip G C~([0, oo), (0, oo)). 
Setting v r (x) := ?; T (M)|fj- It is easily checked that v T satisfies HAGS I (8.3.8)] for all radially 
symmetric test functions <f>. If (f) G C£°(R n x (0,oo)) is a general test function, the function ip, 
defined by 

^ PjT ):=—^—[ (f)(x,T)d?{ n ~ 1 x (3.15) 



Pn(p) 



dB a 



is in C^° r ([0, oo), (0, oo)). By definition of v T and ip, it follows that 

d<p(x, t) 



+ (v T (x), V x 4>(x, r)) u(x)dxdr 



o Jo 



+ ^T^M I "MftMdpdr = 0, (3.16) 



which means that v satisfies HAGSl (8.3.8)] for any test function. Hence, r i-)- u(r) is absolutely 
continuous, and |m'(t)| < ||v r || L2(IR m u) = ||v T ||£2(( 0)OO ) )U ( T )p n ). □ 

Remark 3.5. If follows in particular that, if t K u(r) is absolutely continuous, there exists v 
satisfying (13.111) such that \\v T ||l 2 ((o,oo),«p„) = Kl( r ) for almost all r. By HAFPl Proposition 8.4.5], 
this v T is uniquely determined for almost every r. 

Note the similarity between Lemma [3731 and the theorem cited above: both relate absolute con- 
tinuity and the metric derivative to the solvability and a norm of solutions of a certain weakly 
formulated equation. Both equations will prove useful later. 
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Intuitively, the optimal transport map and the solution of the continuity equation are related: the 
former is displacement of mass, the latter is velocity of mass. This intuition can be made precise 
using HAGSl Propisition 8.4.6], where it is shown that if r (->■ u{t) is absolutely continuous, 

lim^- = v T (3.17) 

h^O h 

for almost every r, where v T is the solution of ( 13.111) such that ||fT||L 2 ((o,oo),M(r)p„) — |w'|(t) for 
almost every r. 

Absolute continuity in X n can now be characterized in terms of absolute continuity in C n and 
L*. That is, 

Lemma 3.6. A curve r i-» (r(r), w(r)) in C n x L x n is absolutely continuous if and only if the curves 
r I—)- r(r) and r i-> w(r) are absolutely continuous in C n and L*, respectively. In this case, 

l(r,u)1(r)=(|r'|^(r) + ^^) 5 (3.18) 

for almost all r. 

Proof. Let r i-> (r(r), u(t)) be a curve in C n x L\. By definition, 

d(r(s),r(t)) < e((r(s),u(s)),(r(t),u(t))), (3.19) 
W 2 (n( S ),«(t)) < y/^g((r(s),u{s)),(r(t),u(t))), (3.20) 

which means that absolute continuity of r i— >■ (r(r), «(r)) implies absolute continuity of r i— >• r(r) 
and r i-)- u(r). 

Conversely, if both r i-)- r(r) and r i— > -u(r) are absolutely continuous, 

Q((r(s),u(s)), (r(f),«(f))) = (d 2 (r( S ), r(t)) + E!M^)T\ 2 

<-\[flrUr^)\(r^ar)i (3.21) 



<2 /Y|rt(r) + ^^Vdr 



which means that r i-> (r(r), u(r)) is absolutely continous. 
Finally, if r i— > (r(r), "u(t)) is absolutely continous, 

|( y|(r) = ^ m Q((r(r + h)Mr + h))Mr)Mr))) 
h-±0 \h\ 



2\ 2 



lim d l r ( T + ; ')' r ( T »j + I ( [im W i(u(t + ft).u(r)) 



r'|^(r) + 



ft 

where, for almost every r, existence of all limits is guaranteed by absolute continuity. □ 
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Together with the results about absolute continuity in C n and L\, this lemma is the main tool to 
study absolutely continuous curves in X n . 



3.2 Constant speed geodesies 

In a metric space S, a constant speed geodesic is by defintion a curve 7 : [0, 1] — > S satisfying 

d( 7 (s),7(*)) = (*-sM7(0),7(l)) (3-23) 

for any < s < t < 1. Note that this is different from geodesies on a Riemannian manifold: 
geodesies are parametrized by constant velocity on the unit interval, instead of parametrized by 
length. By the triangle inequality, it is sufficient to show only 

d(7(a),T(*))<(*-a)d( T (0),7(l)). 

It is also clear that a constant speed geodesic is absolutely continuous, and the metric derivative is 
equal to <i( 7 (0), 7(1)) almost everywhere. The converse is also true: if 7 is absolutely continuous, 
and 

W\ d (t)<d(rf(o),j(i)), 

then 7 is a constant speed geodesic. From this observation, it follows that r >-)■ (r(r), h(t)) is a 
constant speed geodesic if and only if the maps r4r(r) and r 1— > u(r) are. 

Lemma 3.7. A curve r (r(r),M(r)) in C n x L\ is a constant speed geodesic if and only if 

r 1 — ^ r(r) and r 1— > -u(r) are. 

Proof. Suppose first that r 14 r(r) and r h-> m(t) are constant speed geodesies. Then, using Lemma 
13.61 t 1-4 (r(r), u(t)) is absolutely continuous, and 

Kr.-O-IV) = MiW + ^ =d»(r(0),r(l)) + MM^W) ^ 
= £» 2 ((r(0),«(0))?(r(l),«(l))). 

Conversely, suppose that r i-> (r(r), tt(r)) is a constant speed geodesic. Using Lemma [3T6l and 
Jensen's inequality, 

r 1 \u'\ 2 (r) 

g 2 ((r(0),u(0)),(r(l),u(l)))= \ \rf d (r)+ l -^^dr 

(3.25) 



> QVld(r)dr) + ^j\ u '\ W2 {r)dr 

>^ 2 ((r(0),n(0)),(r(l), M (l))) 

which means that all inequalities are in fact equalities. Then r 1— > \r'\(r) and r 1— > \u'\(r) must 
be constant for almost all r. By the second (in)equality, these constants must be d(r(0), r(l)) and 
W / 2(w(0), w(l)), respectively. It follows that r i-> r(r) and r i-> w(r) are constant speed geodesies. 

□ 



19 



With this lemma at hand, it would seem that the constant speed geodesies of d and W2 can be 
studied separately. This is not the case, as X n is a subspace of C n x L\ . After characterizing the 
constant speed geodesies of d and W 2 , the resulting geodesic in C n x L\ will be studied. 

Using the isometry (12.141 ), it is easy to show that the constant speed geodesies of d are given by 

2 

r(r) = ((1 - r)r(O) 1 ^ + rr(l) 1 ^) " +1 > (1 - r)r(O) + rr(l). (3.26) 

It is shown in HAGS I Theorem 7.2.2] that the constant speed geodesies of the Wasserstein metric 
are given by 

«(r)= ((l-r)id + <gj) «(0), (3.27) 

where t^fy is the optimal transport map from u(0) to w(l). This expression has a nice interpretation: 

the transport map tj^l tells where the mass at a certain position has to go to change the profile u(0) 
into u(l). Thus (13.271) is the evolution where all mass travels from its initial position to its destination 
at a constant speed. 

The characterization (12.211) of the optimal transport map, which followes from the radial sym- 
metry, leads to the following characterization of a constant speed geodesic, which is also shown in 
HAGSl (7.2.8)] 

K { l ) = a-r)F^ ) +rF- { l ) (3.28) 

where F u is defined as in (12.201 ). 

As noted above in Lemma [3771 the constant speed geodesies of C n x L\ curves r i-)- (r(r), «(t)) 
with r and u as in (13.261) and (13.271) . By the inequality in (13.261) . (r(r), u(t)) E X" if this is the 
case for r = and r = 1: geodesies between points in X n do not leave X n . Note, however, that 
the geodesies have a peculiar property. Since the inequality in (13.261) is strict for < t < 1 unless 
r(0) = r(l), the support of u{r) will be strictly smaller than B r ( T \, even if w(0) and u(l) are positive 
throughout B r ^ and B r ^, respectively. 

3.3 A-convexity 

As announced above, the concept of A-convexity will be used. A-convexity extends the notion of 
ordinary convexity, which will be equivalent to O-convexity. The number A can be interpreted as a 
measure of how convex a functional is. 

Lacking a linear structure, convexity of a functional on a metric space has to be defined using 
curves. A functional on a metric space (S, d) is said to be A-convex along 7 : [0, 1] — > S if 

<j>(l{r)) < (1 - r)#y(0)) + T0( 7 (1)) - \r{\ - r)d 2 ( 7 (0), 7 (1)) (3.29) 

for all r G [0, 1]. Clearly, it cannot be expected that a functional 4> is A-convex along all curves. 
Inspired by convexity of a function on Euclidean space, which is equivalent to O-convexity along 
straight lines, one usually aks whether a functional is A-convex along geodesies. 

With the results from the previous section in mind, the convexity of $ and q 2 can be studied term 
by term. It will be shown that $ is A-convex along geodesies, and that (r, u) i-)- ^g 2 ((r, u), (s, w)) 
is 1-convex along geodesies. 
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First of all, J 7 is 0-convex along geodesies if the map 



z z n f (z- n ) 

is convex and nonincreasing, as is noted in HAGS I Proposition 9.3.9] and was first shown by McCann 
llMl Proposition 1.2]. Similar to the proof of Lemma |2.5[ this is implied by 

A (/7 = -nz n - l f (z- n ) , (3.30) 

which is an increasing nonpositive function. It will turn out that this property of / also plays a role 
when computing the local slope of T . 

The convexity of W$ can easily be checked in the radially symmetric situation. From (|3.28l) and 
(|2.22l) . it follows immediately that 

W*(u(t),w)= f \{l-r)U \a) + rU^{a)-W-\a)\ 2 da. (3.31) 
■/(o,i) 

A straightforward calculation shows that the map x H- \\x — y\ 2 h 1-convex, as is also shown 
in HAGS I Remark 2.4.4]. Combining this with the above expression for W% implies that w h-> 
^W$(u,w) is 1-convex along geodesies for any u G L\. Note that the special properties of W 2 
for radially symmetric profiles are used in the proof. This is really necessary: in the general case, 
the Wasserstein distance is not 1-convex. A straightforward counterexample is given in [AGS] Ex- 
ample 9.1.5]. Actually, it is shown in HAGS I Theorem 7.3.2] that the opposite inequality holds. 

The convexity of P n can be checked using the isometry (12.141) : A-convexity of P n is equivalent 
to A-convexity of the map 

2n-2 

s H- 9 n s "+ 1 

where 9 n > is a constant depending on the dimension. Note that the exponent is always between 
and 2. The second derivative of this map is 

_ (2n-2)(n-3) ^ 
(n + \y 

Using HAGSi Remark 2.4.4], this means that P n is 0-convex if n > 3. Note that P 2 is not A-convex 
along geodesies for any A since the second derivative is not bounded from below. However, if s is 
bounded away from zero, the second derivative is bounded from below, and P n is A-convex, with A 
equal to the infimum of the second derivative. As was shown in the previous section, the radius r is 
bounded away from in any sublevel of which means that restricting the problem to a sublevel 
means that P n is A-convex for some A < 0. Note that a similar argument can be used for n > 3 to 
obtain A > on sublevels. 

Finally, again using the isometry (12.141 ) and RAGS I Remark 2.4.4], s i-> |d 2 (r, s) is 1-convex 
along geodesies for all r E C n . 

Collecting the above results, one finds 

Lemma 3.8. The map (s, w) i->- |f? 2 ((r, u), (s, w)) is 1-convex along geodesies. Moreover, if n > 3, 
$ is 0-convex along geodesies. In case n = 2, $ is A-convex along geodesies with A < on its 
sublevels. 
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In particular, Assumption B.ll holds. That is, given (s, w), (r(0), u(0)), (r(l), u(l)), the map 

(r,w) ^$(r, W ) + ^-f ? 2 ((s,w;),(r,n)) (3.32) 

is convex along the constant speed geodesic r i-)- (r(r), u(r)), defined by (13.261) . (13.271) . 

One can ask if the lack of convexity in case n = 2 is due to the odd shape of the geodesies 
of d. Obviously, the convexity of P n is much better along linear interpolants than along geodesies: 
even if n = 2, it is immediately clear that P n is 0-convex without having to restrict to a sublevel of 
$. Unfortunately, d 2 is not 1-convex along linear interpolants. Moreover, the profit from replacing 
geodesies with linear interpolants is somewhat disappointing: even in case n > 3, it is not possible 
to obtain A-convexity for A > 0, which would be the metric equivalent of strict convexity. Therefore, 
it does not seem to be beneficial to use other curves instead of geodesies. 

Having shown that Assumption 13.11 holds. [AGS, Theorems 2.4.15 and 4.0.4] can be applied. 
Since the set of (r, it) with $(r, u) < +oo is dense in X n , no condition on the initial value is needed 
anymore. 

Theorem 3.9. For any (r , u ) e X n , 

• There exists a unique mimimizing movement £ !->■ (r(£), u(t)) for $ starting at (r , u ). 

• (r(t),u(t)) is also the unique generalized minimizing movement for $ starting at (r , uq). 

• (r, it) is a locally Lipschitz curve of maximal slope with \d$\(r(t), u(t)) < +oo for t > 0. 

• (r, u) is the unique solution of the evolution variational inequality 

~Q 2 ({r{t),u{t)), (s, w)) + ^g 2 ((r(t),u(t)), (s, w)) + ®(r(t),u(t)) < w) (3.33) 

for all (s, in) G X n such that $(s, w) < +oo and almost all t > 0. 

• For (r , Uo) and (so, Wo), the minimizing movements (r, u) and (s, w) for $ starting at (r , ito) 
and (so, wo), respectively, satisfy 

e ((r(t),u(t)),(s(t),«;(t))) < e- A ^((r , Mo ),(s ,wo)). (3.34) 
If $(r , u ) < +oo, the right metric derivative 

|(r,<|( t): =lim g « r (' + / ')'"( f t'' )) ' (r( ' ) '" ( ' ))) . (3-35) 

and the equation 

d -$(r(t), M (t)) = -|a$| 2 (r(t), M (t)) = -|(r,<| 2 (t) = -|fl*|(r(t), «(t))|(r, (3.36) 



holds for all t > 0. If, additionally, |<9$|(r , u ) < +oo, (13.361) also holds for £ = 0. 

Remark 3.10. The definitions of curve of maximal slope and |<9$| will be presented below. It will 
turn out that (13.361 ) is stronger that the statement that (r(£), it(£)) is a curve of maximal slope 

Unfortunately, as A < 0, the contraction property does not give a lot of information about the 
asymptotic behaviour of the minimizing movement. 
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4 Curves of maximal slope 



One of the conclusions of Theorem |3.9| is that for any (r , u ) e X n , a curve of maximal slope for $ 
starting from (r , «o) exists. In this section, the concepts of local slope, denoted by |<9$| and curve 
of maximal slope will be introduced. Additionally, the local slope of $, which may be infinite, will 
be computed in this section. 

A curve of maximal slope can be regarded as being the metric equivalent of a gradient flow. As 
before, the inspiration is the situation in R™. The equation x(t) = —V<p(x(t)) flow for the gradient 
flow of some smooth : R n — > R is equivalent to 

^(x(t)) = -i|V0(x(t))|-i|i(t)|. (4.1) 

Using Young's inequality, 

±<Kv(t)) > -\V<P(y(t))\\m\ > -lW(y(t))\ 2 - \\m? (4-2) 
for any smooth curve y. Therefore, a gradient flow in R™ can be characterized by the inequality 

A0( x(t)) <_I|V0(x(t))| 2 -i|x(t)| 2 . (4.3) 

If it is possible to redefine | V0| and |x| in a metric space such that the chain rule for <p o x still 
holds, this inequality can be used to define a gradient flow in a metric space. A generalization for 
\x\ has already been found: it is the metric derivative defined in the previous section. A possible 
generalization for | V0| is the local slope. The local slope of a functional <p : S — > R is defined by 

M(»0=limBup (4.4) 

w ^ v d(v,w) 

Note that the local slope is in some sense one-sided: it only measures how fast the value of <p 
decreases near a point. This makes sense, since the functionals that are considered are usually only 
lower semicontinuous. See [AGS, §1] for more details. 

The local slope of $ can be computed using methods from the proof of HAGSl Theorem 10.4.6]. 
The main ingredient is HAGSl Lemma 10.4.4], which describes the behavior of T{u) as u is pushed 
forward by a sufficiently smooth map r. 

Lemma 4.1. Let u e L\, r e L 2 (R n , it) be a radial vector field, f > and suppose that 

1. r is differentiable w-almost everywhere, and r t := (1 — r)r + rid is w£ n -injective with 
| det Vr T | > it£ n -almost everywhere for any r 6 [0, r], 

2. || V(r - id)!!^^^) < +oo 

3. There exists a constant C > such that V^, z 2 : f(zi + z 2 ) < C (1 + f(zi) + f(z 2 )) 

4. J r ((r T ) # M) < +oo. 
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Then the map r n> r 1 (J r ((r T ) # w) — T(u)) is nondecreasing in [0, r], and 

+ ex) > lim Jr((lv ) #M ) ~ = _ f /( u (x))trV(r(a;)-x)dx (4.5) 

It is easy to check that 2; 1 — >■ 2; log z satisfies this satisfies Condition ©, also known as the 'dou- 
bling' condition. Moreover, it does not limit the growth of / too much. For instance, the functional 
2 H- 2 2 also satisfies Condition ©. The conditions on r are a bit technical, but there are two im- 
portant examples of vector fields that satisfy the requirements: smooth, injective vector fields and 
optimal transport maps between profiles with copmact support. The latter follows from HAGS1 The- 
orem 6.2.7] about regularity of optimal transport maps, and w£ n -essential injectivity follows as 
in the beginning of the proof of HAGSl Proposition 9.3.9]. Note that the apprximate differential in 
HAGSl Theorem 6.2.7] is not needed if r is a transport map between compactly supported probability 
measures. 

Before the main theorem can be proven, a refined version of HAGSl Lemma 10.4.5] is needed. 
This lemma shows that, if / is smooth enough, a weak integration by parts formula holds. 

Lemma 4.2. Let (r, u) E X n , and r and f be as in Lemma 147X1 In addition, suppose that 

1. f{u)\ Br eW^{B r ), 

2. r is of bounded variation. 
Then 

f(u(x)) div(r(x) — x) dx < — / Vf(u(x)) ■ (r(x) — x) dx + (s — r) f{u{r))P n {r) (4.6) 



13 r 



Proof. Without loss of generality, it can be assumed that r = 1. As in the proof of [AGS, Lemma 
10.4.5], the main ingredient is an estimate for the distributional divergence of a BV vector field -q. 
However, now using HAFPl 3.92], the boundary term at dB r should also be taken into account. If 
w E C^°(IR n ) is nonnegative and radially symmetric, its trace on dB r is constant, say w(x) — w if 
|x| = r. Then 

/ w(x) div(rj(x)) dx < — / Vw(x) • rj(x) dx +r)(r)wP n (r) (4.7) 

J Bf J Br 

if the distributional divergence of rj is a nonnegative measure. For bounded 77, the same inequality 
holds for nonnegative and radially symmetric w such that w\b t E W l,1 (B r ) and w = outside B r . 
Applying this for w = f(u), rj := r, 



f(u(x)) div(r(x)) dx < — Vf(u(x)) ■ r(x) dx + r(r)f(u(r))P n (r) 

B r J B r 



(4.8) 



< - / V/(«(x)) ■ r(x) dx +sf(u(r))P n (r) 



Standard integration by parts yields 



f(u(x))div(x)dx = — / V/(u(x)) • xdx+rf(u(r))P n (r), (4.9) 

B r J B T 

which, together with (14.81) , yields the result. □ 
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Comparing this result to HAGSl Lemma 10.4.5], the approximating sequence is left out. 
This can be done because the lemma will be applied to an optimal transport map between radi- 
ally symmetric measures. From (|2.21l) . it can be seen that t™ is always monotone. Hence, the map 
t™(x) := £^(a?)o, f° r which the lemma will be applied, is of bounded variation. 

Apart from radial symmetry, the only difference between this proof and the proof of HAGSl 
Lemma 10.4.5] is from the boundary term in (14.71) . In the following computations, more boundary 
terms will appear in a similar way. Loosely speaking, VF 1,1 (M n ) functions in HAGSl §10.4] are 
replaced by functions that are W 1 ' 1 when restricted to B r . In the language of BV functions, this 
means that the distributional derivative consists of an absolutely continuous part and a jump on 
dB r . Obviously, it is not important that the domain is a ball when performing these calculations: the 
same calculations can be done on much more general domains. 

Since the calculation is basically integration by parts, it is not surprising that boundary terms 
appear when restricting to a bounded domain. Through the calculations that follow, however, the 
boundary term will give the osmotic term that was discussed in the introduction. This will become 
clear in the proof of the following theorem that gives the local slope \d$\ of 

Theorem 4.3. Let (r, u) 6 X n be given. Then $ has finite local slope at (r, u) if and only if 

f(u)e ^((Cr), P n ) with 

V/(«) 



e L 2 m n ,uP n ). 



u 



(4.10) 



In this case, 



\d$\(r,u) 



n — 1 



/Mr)) 



PJr) + k 



V/(u) 



u 



(4.11) 



L 2 (R",u) 



Proof. Suppose first that |d$|(r, u) is finite. Let £ be a radial vector field, and define £ by £(ar) 
£(|x|) A. Similar to the proof of HAGSl 10.4.6], note that 



W 2 (u, (X T )#u) < r||£|| L 2 (IR n iU) = r||£|| L 2 ((0iOo)iUPn) 



(4.12) 



where X T (x) = x + t£(x). Using the definitions of \d&\ and metric derivative, followed by this 
inequality and (13.51) . 



f(u(x)) div£(x) dx 



^ Aim mX T (r),(X T )#u)-<S>(r,u)\ | |P w (X r (r)) - P w (r)| 

~ tIO t t 



< |3$|(r, u) \{Xr{r), (X T ) # u)\' (0) + |£(r)|i^(r) 



< \d$\{r,u) |£(r)| 2 P n (r) + 



ll?l| 2 



L 2 ( (0,oo),-uP n ) 
K 



(4.13) 



whenever £ is smooth enough to apply Lemma |4~T1 with r = £ + id. For £ = — £, this yields 

nl /(»(./ ))d.r < \d$\(r,u) I r 2 P n {r)+ H ld ^(( '°°)^) » +rP ^ r ) j (4 .i 4) 



K 
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which means that f(u) is integrable. Next, apply (14.131) for smooth, compactly supported radial £ 
to obtain 



f(u(x)) div^(x) dx 



< \d$\{r,u) |e(r)| 2 P„(r) + 



lien 2 



L 2 ((0,oo),«P n ) 



< [|9$|( r , u )(i + p n ) + J*(r))||£| 



C°((0,oo))- 



(4.15) 



Noting that this inequality only needs to be checked for radial £, Riesz' theorem HAFPl 1.54] implies 
that x i — y f(u(x)) is of bounded variation. Hence, 



(4.16) 



f(u(x)) div dx = - I dDf(u(x)) . 

Restricting to £ with £(r) = 0, 



£dD/(«) 



f(u(x)) div £(x) dx 



JJ r 



< ia$|(r,u) 2 



L 2 (K",m) 



(4.17) 



Using duality, this means that Df(u) is absolutely continuous on i3 r . In particular, /(w)|s r G 
PK 1,:L (-B r ). Since u = outside _B r , it follows that the singular part of Df(u) must be concen- 
trated on dB r . By radial symmetry, it follows from flAFPl Proposition 3.92] that Df(u) only consist 
of an absolutely continuous part concentrated on B r and possibly a jump part accross dB r . More 
precisely, 



f(u(x)) div £(x)r dx 



13 r 



x 



V/(u(x)) • £(x) dx + / f(u(x))£(x) ■ - dH n -\x) 



8B r 



dp 



f{u{p))t{p)P n {p) dp + f{u{r))ti{r)Pn{r). 



(4.18) 



Using Lemma S?U and (14.181) 



lim «XAr)AXr)**)-*lr,u) = 

riO T 



f(u(x)) div £(x) dx 



P' n {r) - f(u(r))P n (r) ) £(r) + / V/(«(x)) • £(x) dx . 

(4.19) 



Estimating the left hand side as in (|4.13l) . 



PJr) 



f(u(r)) P„(r)£(r) + 



Vf(u(x)) 



B r U[X) 



£(x)u(x) dx 



< \d$\(r,u) |£(r)| 2 P n (r) + 



|£l|2 



L 2 (R", 
K 
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Finally, using another duality argument and evaluating 



P n (r)> 



n — 1 



f(u(r)) 



Br 



Vf(u(x)) 



U[X) 



u(x) dx < |<9$|(r, u). 



(4.20) 



Conversely, assume f(u) | b t £ W 1 ' 1 (B r ) 



V/(u) 



e L 



(4.21) 



and let (s, w) G X n such that $(s, w) < $(r, tt). As noted above, the map £™ (x) = t™(x) A, wh ich 
is the optimal transport map from m to in R™, is sufficiently regular to apply Lemmata 14. II and 
UThen 



$(s,w)-$(r,u)>(P n (.s)-P n (r))- / /(«(*)) trV(t£(ar) -a) da; 



> 



/(u(r)) (a - r) P n (r) + / V/(u(z)) • (t™(x) - x) dx 



> - 



n — 1 



f(u(r)) 



Pn(r) + K 



v/(«) 



L 2 (R™,«) 



|s - r| 2 P n (r) + 



Wi(u,w) 



(4.22) 



using convexity of r i-> P n (r). The inequality in (14.1 II) left to prove now follows by writing a Taylor 
expansion for the integral from (12.1 II) . and taking the limit superior for (s, if) — )■ (r, u). □ 



Note again that boundary terms appear, most importantly in (14.181) . As explained above, this is a 
direct consequence of restricting u to a bounded domain. The fact that the boundary of the domain 
moves is not important here yet. In (|4.17l) . it becomes clear that the boundary term disappears when 
the perturbation is such that the mass at the boundary does not move. Studying the calculations a bit 
closer, the boundary term can be 'estimated away' if £(r) < 0. Hence, the boundary term still plays 
a role when the boundary is fixed, but disappears from the final outcome. From a modelling point 
of view this makes sense: if the cell membrane would be fixed for some reason, the osmotic force 
would still be there. The only thing that has changed is that te membrane does not react to forces 
anymore. In a more general setting, a fixed boundary can be regarded as a moving boundary with 
infinite resistance to force: it can be seen from the calculation that all forces are there, but disappear 
because the boundary is unable to react to them. 

The section is concluded with a lemma that essentially characterizes the minimal Frechet subd- 
ifferential of $. 
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Lemma 4.4. Suppose that (r(r),-u(r)) is an absolutely continuous curve in X n , and v(t) is the 
solution of (|3.1 II) with ||^ t ||l 2 ((o,oo),«(t)p„) — l w 'l( r ) f° r almost every r. For any r such that 
|(9$| (r(r),«(r)) < +00, r i-)- $(r(r),M(r)) is differentiable, and (13.171) holds, 



U 1 - /Mr(r),r))j Pn(r(r))r'(r) + jf^ v T (p)P n (p) dp. 

(4.23) 



d Mr(r)Mr)) y 



Proof. For /i > 0, let be the optimal transport map from u t to u t+h . Then, as in the proof of the 
Theorem 14 .3 [ 



$(r(r + /i), u(r + h)) - $(r(r), w(r)) 

n — 1 



> 



r(r) 



- /(«(r(r),r)) (r(r + h) - r(r)) P„(r(r)) 



+ 



r ^df(u(p,T)) <T _ 



(tJ T ; h) (p)-p)P n (p)dp (4.24) 



i(T) 



'0 #P 

Using (|3.17l) . dividing by /i, and taking limits for /i 4 yields 

A$( r(r )Xr)) > - /(u(r(T),r))) P n {r{r))r\r) + jf^ ^ ( ^' T)) v T {p)P n {p) dp 

(4.25) 

where f is characterized by Remark [331 The converse inequality follows by studying the left deriva- 
tive. □ 



5 Weak Solutions 



A natural question is whether a maximal slope curve obtained in Theorem 13 .9| is a solution of (11.71 ), 
and in what sense. Bearing in mind Lemma |44l it is expected that the argument from the beginning 
of the previous section justifying the definition of a curve of maximal slope can in this case be 
reversed to obtain a weak form of (11.71) . Of course, the choice for / should be kept in mind: 



u t = «A/(u), 
= u{r)r'{t) 

77,-1 



r' 



+ f(u(r)) 



for x G B r (t), t > 0, 
for t > 

for t > 0. 



(5.1) 



It is possible to formulate this problem in a neater way. Note that 

V/(u) = V (uf'(u)) - V/(«) = f(u)Vu + uVf\u) - f{u)Vu = uVf(u) (5.2) 
if / is sufficiently smooth. Therefore, the first equation in (15.31) can be replaced by u t = kA/(m). 
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Similarly, the boundary condition for u can be rewritten: 

u t = ftdiv (uV/'(u)) , 
= r'(t), 



p=r 



n — 1 



+ /(«(r)) 



for x G -B r (t), £ > 0, 
for t > 

for £ > 0. 



(53) 



This way of writing the problem makes the structure clearer: particles have velocity equal to V/'(u), 
and the velocity of the particles at the boundary should match the normal velocity of the boundary. 

In order to show that a curve of maximal slope can be regarded a weak solution of (15.31) . some 
integral identities are derived that characterize solutions when smoothness is assumed. By smooth, 
it is meant that t !->■ r(t) is continuously differentiable, and (x,t) !->■ u(x,t) is twice continuously 
differentiable on the domain 

{(p,t) :0<p<r(t)-t>0} 

Moreover, it is shown that $ is a Lyapunov functional for (15.31) . As before, when applying the 
divergence rule, u(t) is interpreted as a function on W 1 . 

If u is smooth in the sense explained above, and ip G C£°(R n x (0, oo)) is radially symmetric, 



J B 



u t (x, t)ip(x, t) dx dt 



r(t) 



u t (x, t)(ft(x, t) dx dt 



o J B r 

+ f u(r(t),t)<p(ty(t)P n (r(t))dt, 



oo 



(5.4) 



where 0(t) is the value of </?(., t) on dB r u). Moreover, 



o J B 



u t (x, t)ip(x, t) dx dt = —k 



r(t) 



B r(t) 



Af(u(x, t))ip(x, t) dx dt 
u(x, t)Vf(u(x, t)) ■ Vip(x, t) dx dt (5.5) 
0(t)P n (r(t)) dt 



dp 



p=r 



Combining these two identities with the Neumann boundary condition yields 



J B r 



u(x, t)(p t (x, t) dx dt 



(*) 



o Jb 



Vf(u(x,t)) ■ V<p(x,t)dxdt. (5.6) 



r(t) 



By a similar computation, for any radially symmetric ip G C^°(R n x (0, oo)), 

POD 

ipt(x,t)dxdt = - / 4>(t)r'(t)P n (r(t))dt 



•* Br(t) 
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- f{u(r{t),t))) $(t)P n (r(t))dt, 



(5.7) 



where ip is defined analogous to (p. From these identities, it follows that $ is indeed a Lyapunov 
functional for (15.31) : 



dt 



Pn(r(t))+ / f(u(x,t))dx 



K(r(t))r'(t)+ / f'(u(x,t))u t (x,t)dx + f(u(r(t),t))r'(t)P n (r(t)) 
K ' [r(l)] f(u(r(t),t)))r'(t)P n (r(t)) 



Pn(r(t)) 

+ K f'(u(x, t)) div (u(x, t)Vf'(u(x, t))) dx 



r(t) 



+ u(r(t),t)f(u(r(t),t)y(t)P n (r(t)) 



1 \ 2 

n — 1 



/(u(r(t),*)) P n (r(t)) + K \Vf'(u(x,t)\ 2 u(x } t)dx 
+ /'(w(r(t), t))P„(r(*)) ( u(r(t), t)r'(t) + k -ff(u(p, t)) 

TW - /(M(r(tM)) ) P " m) - K L 



since 



u [/' o u]' 



(5.8) 



If the integral identities (15.61) . (15.71) . together with monotonicity of $ along a trajectory are ac- 
cepted as a definition of weak solution of (15.31) . the following theorem simply states that curves 
of maximal slope of $ are weak solutions of (15.31) . Loosely speaking, this means that (15.31) can be 
formulated as a gradient flow. 

Theorem 5.1. Suppose that t (-> (r(t),u(t)) E X n is continuous. Then (r(t),u(t)) is a curve of 
maximal slope for $ if and only if $(r(t) , w(t)) is decreasing, and (|5.6I) and (15.71) hold for all radially 
symmetric test functions 4>,ip E C%°(M. n x (0, oo)). 

Proof. Assume first that (r(t),u(t)) is a curve of maximal slope for $. By definition, this means 
that <p{t) : = $(r(t), u(t)) is decreasing, and t h-» u(t) and £ i-> r(t) are both absolutely continuous. 
Then there exists a function v satisfying (13.1 II) and ||i> r ||L 2 ((o,oo),uP n ) = \ u '\(t) f° r almost every 
t > 0. Hence, as in the proof of Theorem 13 .41 v T (x) := v T (\x\)-^ solves 

rOC I* rOO r 

I I (f t (x,t)u(x,t) dx dt = — / / V<^t • v(x, t)u(x, t) dxdt (5.9) 

Jo JB r{t) Jo JB <t) 

for any radially symmetric ip E C^(W l x (0, oo)). Using Lemma |4~4l 

= ( — - f(ut(r t ))] P n (n)^ + T (T) 9 ^ P,T)) v T (p)P n (p) dp (5.10) 



n J dt J dp 

30 



for almost all t > 0, where (f>(t) := $(r(t), u(t)). On the other hand, 



m < 



\(r,u)'\ 2 (t) \d$\ 2 (r(t),u(t)) 



\r'(t)\ 2 P n (r(t)) + 



n — 1 



r(t) 



ll V *ll^(R", M (f)) ^ 



/(«(»■(*),*)) 


Pn 




Vf(u(x,t)) 


2 


R n ,u(t))) 


u(x, t) 


L 2 (l 





(5.11) 



for almost all t > 0. By Young's equality, it follows that for all £ > such that both hold, 

n — 1 



r'(t) 



r(t) 



V t = -AC 



+ f(u(r(t),t)), 

«(t) 



Substituting the latter in (f5T9b yields (I5T61) . Together with (|377T) . the first implies (15771) . 

Conversely, note that (|5.6I) is equivalent to (13.1 II) with a given expression for v(t) substituted. 
Therefore, t !->■ is absolutely continuous with respect to the Wasserstein metric. Similarly, by 
Lemma 1531 1 \-t r(t) is absolutely continous, and gives and expression for r'(t) almost everywhere. 
Then Lemma l4~4l yields 



n — 1 



f(u(r(t),t)))P n (r(t)y(t) + 



^ df(u( P ,T)) 

dp 



v T (p)P n (p)dp (5.12) 



r(t) 

for almost every t > 0. 

Having explicit expressions for ^ and ^ at hand, it follows from Young's inequality that t i— >■ 
(r(t), is a curve of maximal slope. □ 



Note that the splitting of the problem from the introduction plays an important role: the equa- 
tions related to absolute continuity give the equations relating the evolution of u and r to different 
velocities, whereas the gradients give the velocity in terms of the current state. As argued above, the 
divergence formulation (15.31) . also has this structure. 



6 Varying permeability 

As stated in the introduction, two variants of the cell swelling model could be considered. Apart 
from the above, which corresponds to a fixed permeability of the membrane, an alternative model 
based on aquaporins could be studied. In this section, the model will be adjusted to reflect this, and 
it will be shown that the gradient flow approach above also works for the resulting problem, with 
only slightly different results. 
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6.1 Adjusting the model 



If it is assumed that the number of aquaporins on the cell membrane is fixed, and the permeance of 
an aquaporin is also fixed, the total permeance of the membrane is fixed. By radial symmetry, the 
permeability of the membrane is equal to p"^)) > U P t0 a multiplicative constant. This means that 
the new equation for the normal velocity of the membrane will be 

As before, all but one of the parameters can be made equal to 1 . The one remaining parameter will 
again be k. In this case, however, the scaling factor for t also depends on r(0). 

From a modeling point of view, there is only one reasonable way to adapt the gradient flow 
approach to this new situation: since permeability basically determines how hard it is to move the 
membrane, the metric for the r-coordinate should be modified. Remembering that d arose from 
minimizing 



\\Vt\\L*{dB r{t) )&t, 

it seems reasonable to integrate some other norm of v t . It turns out that choosing the reciprocal of 
the permeability, that is, P n as a weight in the L 2 -norm of v t is the right choice. Therefore, define 

d(r ,n) := ! P n {p) dp = C n (B ri AB ro ) = u n |< - r "| . (6.2) 
Obviously, the metric for C n x L\ will be the same combination of d and W 2 as before: 



g((r,u),(s,w)):= ( d 2 (r, ,) + ^ w) ) . (6.3) 



6.2 Implications of changing the metric 

Note that the new definition does not change the topology on C n , which means that most of the 
results still hold. Some of the finer properties, however, have changed. In this section, the results 
from the previous sections will be adapted to reflect the new model. 

First of all, since the metric has changed, the isometry from C n to [0, +oo) has changed as well. 
The isometry now reads 

t n : C n ->• [0, +oo) : r \-> u n r n . (6.4) 

All other results from Section [2] are the same. 

The metric derivative |/|^ is of course also different. Absolute continuity in [0, +oo) and C n are 
of course still equivalent, but 

KU(r) = P„(r(r)) 

for an absolutely continuous curve r !->■ r(r). The obvious adaptation of Lemma [33] also holds. 



dr(r) 



dr 



(6.5) 



32 



Lemma 6.1. Let r i— >• r(r) be a curve in C n . Then r(r) is absolutely continuous if and only if it is 
continuous and there exists a function g such that g-(r)P n (r(r)) G L^^O, oo)), 



r)^(r(r),r)P n (r(r))dT 



(6.6) 



for all ip G C£°([0, oo) x (0, oo)). In this case, |r'|^(r) = | g (r)|P n (r(r)) for almost every r > 0. 
As before, the isometry t n can be used to characterize the geodesies of d: 

r(r) = ((1 - r)r(0)" + rr(l) n )» > (1 - r)r(0) + rr(l), (6.7) 

which means that combining geodesies for r and u will give a geodesic for g. 

A more substantial difference is the change in A-convexity of P n . By the same argument as 
before, A-convexity of P n is now equivalent to A-convexity of the map 



s i — y 9„r 



n-l 



n. 



(6.8) 



where 6 n > is a constant. Before, P n was 0-convex for n > 3, but this is not the case anymore: 
for every n, the situation is similar to the situation n = 2 before: P n is not A-convex globally, but 
only A-convex with A < on sublevels of $. 

Together with the observation that s i-> |d 2 (r, s) is still 1-convex for all r, HAGS I Theorems 
2.4.15 and 4.0.4] still apply. Since A is still negative, no stronger contraction results follow from 
these theorems. 

Of course, the local slope <9$ with respect to g should be different. Since the proof of Theorem 
14 .31 largely deals with the mass component, it only needs a small modification. 

Theorem 6.2. Let (r, u) G X n be given. Then $ has finite local slope at (r, u) if and only if 

f{u) G W l > l {B r ) with 

V/(«) 



eL 1 



u 



In this case, 



|9$|(r, u) 



n 



1 



/(«(r)) 



K 



v/(«) 



u 



(6.9) 



(6.10) 



L 2 (R",u) / 

Proof. The proof is similar to the proof of Theorem 14.3 1 The estimate (|4.13l) now should be 



f(u(x)) div^(x) dx 



< \d$\(r,u) 



r)?P 2 Ar) + 



L 2 ((0,oo),uP„ ) 
K 



+ Wr)\P'(r) 



Proceeding as in the proof of Theorem 14 . 3 1 no w yields 



(6.11) 



PJr) 



f(u(r)) P n (r)iP(r) + 



Vf(u(x)) 



Br- 



u(x) 



■ £(x)u(x) dx 



< \d®\(r,u) 



r)\ 2 Pl(r) + 



Il 2 (R", 

K 
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which, with a duality argument, gives one of the implications. 

The converse implication is also shown as in the proof of Theorem 14.31 but with the Cauchy 
inequality applied differently: 



w) - $(r, u) > {P n {s) - P n {r)) 
'PLir. 



> 



> 



n — 1 



f(u(x))div(t™(x)-x) dx 
f(u(r)))(s-r)P n (r)+ [ Vf{u(x)) ■ (t™(x) - x) dx 

v/(«) " 



r 



- f(u(r)) 



K 



U 



L 2 {W"-,u) 



Wi(u,w) 



Is' 



which, after taking limits, again concludes the proof. 



(6.12) 
□ 



With expressions for the metric derivative and local slope at hand, the proof of Theorem 1511 can 
be repeated. In this new setting it will give 



P n {r{t))r'{t) 



n — 1 



for almost all t > 0, which means that (|5.7I) becomes 



J B 



it) 



dip(x, t) 
dt 



dx dt 



n — 1 



f(u{r(t),t)) (6.13) 



/(«(r(t),t)) U(r(t),t)dt. (6.14) 



The equation for u is, as expected, not changed. 
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